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ABSTRACT 

Newton's  method  for  generalized  equations  (Josephy  [3])  has  been  applied  to 

A 

the  economic  equilibrium  problem  of  the  Project  Independence  Evaluation  System 
(PIES)  Energy  Model, (Josephy- The  resulting  algorithm  involves  solving  a 
sequence  of  linear  complementarity  problems.  Lemke's  complementary  pivot  algo- 
rithm  is  used  for  this  purpose.  In  this  paper , -we-shew, that  the  linear  comole- 
mentarity  problems  will  be  copositive  plus  when  the  negative  of  the  elasticity 
matrix,  -e,  of  the  consumer's  quantity  vs.  price  relation  has  the  following 
properties:  (J.)  positive  diagonals,  (£^)  negative  off-diagonals,  and  (444) 

strict  diagonal  dominance.  These  conditions  are  satisfied  for  Hogan's  example, 
(Hogan — Thus,  Lemke's  algorithm  will  either  converge  to  a  solution  or  show 


that  no  solution  exists.  Under  the  conditions  of* Josephy 


J 


Theorem  1  tJkU 


S' 


solution  to  the  linear  complementarity  problems  will  always  exist.  Hence, 

Lemke's  algorithm  can  be  used  when  the  conditions  of  the  Theorem  1  of  Josephv 

-f^T^are  satisfied.  — " 

- 
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Significance  and  Explanation 

Newton's  method  for  generalized  equations  (Josephy  [3])  has  been 
used  to  compute  solutions  to  an  economic  equilibrium  problem  found  in 
Hogan  [2]  (Josephy  (4]).  This  economic  equilibrium  problem  is  a  simpli¬ 
fied  version  of  the  Project  Independence  Evaluation  System  (PIES) 

Energy  Model  developed  by  the  Federal  Energy  Administration.  The  compu¬ 
tational  algorithm  used  at  each  iteration  of  Newton's  method  in  Josephy 
[4)  is  Lemke's  complementary  pivot  algorithm.  In  this  paper,  we  give 
the  analytical  justification  for  the  success  of  Lemke's  algorithm  in 
the  computation  of  the  economic  equilibrium  for  Hogan's  problem. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


HOGAN'S  PIES  EXAMPLE  AND  LEMKE'S  ALGORITHM 


Norman  H.  Josephv 

1.  INTRODUCTION 

Newton's  method  for  generalized  equations  (Josephy  [3])  has 
been  applied  to  the  problem  of  computing  the  economic  equilibrium 
of  the  Project  Independence  Evaluation  System  (PIES)  Energy  Model 
(Josephy  [4]).  In  this  paper,  we  show  that  the  linear  complemen¬ 
tarity  problems  arising  when  Newton's  method  for  generalized 
equations  is  applied  to  Hogan's  PIES  example  (Hogan  [2])  are 
copositive  plus. 


We  recall  in  section  2  the  structure  of  a  competitive  market 


equilibrium.  Hogan's  example  is  described  in  section  3.  The 


application  of  Newton's  method  for  generalized  equation  to  Hogan's 


example  is  given  in  section  4,  where  the  resulting  linear  comple 


mentarity  problems  are  shown  to  be  copositive  plus 


2 .  EQUILIBRIUM  IN  A  MARKET  ECONOMY 


The  competitive  market  which  underlies  the  PIES  model  consists 


of  two  classes  of  agents,  suppliers  and  consumers,  and  two  classes 


of  goods,  factors  of  production  and  consumable  goods.  The  suppliers 
faced  with  a  perceived  demand  for  consumable  goods,  convert  factors 


of  production  into  consumable  goods,  sustaining  a  cost  for  conversion 


to  and  delivery  of  consumable  goods,  and  charging  a  price  for  those 
goods.  The  consumers  purchase  consumable  goods  at  levels  dependent 
upon  the  prices  of  all  goods. 
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Suppose  there  are  n  consumable  goods.  Let  q  e  and  p  e 
where  is  the  consumption  level  of  the  ith  good,  and  p^  is  the  price 
of  the  ith  good.  For  a  perceived  demand  q  of  consumable  goods,  the 
supplier  will  charge  prices  no  lower  than  p  (q)  e  Rn  .  The  consumer, 

O  + 

desiring  level  q  of  consumable  goods,  is  willing  to  pay  no  more  than 
prices  pD(q)  e  R^  .  A  market  equilibrium  is  a  quantity  vector  q  and 
a  price  vector  p  such  that  p  =  PQ(q)  =  ps (q)  .  That  is,  at  price  p, 
the  market  for  all  goods  will  clear.  The  producers  will  supply  at 
price  p  a  level  of  consumable  goods  q  which  the  consumers  are  willing 
to  purchase  at  price  p.  The  supply  relation  pg  is  typically  determined 
by  the  solution  of  a  cost  minimization  problem  modeling  the  engineering/ 
technological  processes  involved  in  conversion  of  factors  of  production 
to  delivered  consumable  goods.  The  demand  relation  pQ  is  tradition¬ 
ally  a  behavioral  model  econometrical ly  determined  from  historical 
data. 


For  the  case  of  n  consumable  goods,  the  equilibrium  quantity 
q*  and  price  p*  and  factor  levels  x*  for  the  competitive  market  with 
supply  modeled  as  a  linear  program 
minimize  <c,x> 


Ax  =  q 
Bx  =  b 


x  =  0 


total  cost  of  production 
demand  requirements 
non-demand  constraints 
non-negative  factor  levels 


and  a  log- linear  consumer  demand  model 

„  n 


in(q  /q.)  =  l  e  .  tn(p  /p  )  ,  i=l,...,n 
11  j«i  U  i  J 
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where  e  is  the  elasticity  matrix  and  q^  ,  are  fixed  reference 

values  of  ,  that  is,  pD(q°)  =  p°  ,  is  a  triple  (q*,p*,x*)  >  0 

satisfying  the  following  equilibrium  conditions. 


(EQ.  1) 


(EQ.  2) 


(EQ.  3) 


x*  solves:  minimize  <c,x>  subject  to 
Ax  =  q*  ,  Bx  =  b  ,  x  =  0  . 

For  some  multiplier  s*.  (p*,s*)  solve  the  dual 
problem:  maximize  <q*,p>  ♦  <b,s>  subject  to 
ATp  ♦  BTs  =  c  ,  (p,s)  =  0  . 

0  n  q 

tn(q*/qV)  ■  I  e. .  in(p*/p.)  ,  i=l,...,n  . 

1  1  J  J 


3.  THE  HOGAN  PIES  EXAMPLE 


As  an  illustration  of  PIES,  Hogan  [2]  described  a  simplified 
situation  which  included  the  major  aspects  of  the  PIES  model.  Figure  1 
illustrates  the  structure  of  Hogan's  example.  The  factors  of  production 
are  coal,  crude  oil,  steel  and  capital.  The  consumable  goods  are  coal, 
light  oil  and  heavy  oil.  There  are  two  coal  mining  regions,  two  crude 
oil  drilling  regions,  two  oil  refineries,  and  two  demand  regions.  Coal 
can  be  mined  in  region  i  at  level  c  ,  j  =  1,3,  at  three  different 

levels  j  for  differing  production  costs.  Oil  can  be  drilled  in  region 
i  at  level  o^.,  j  =  1,2,  at  two  different  levels  j  for  differing 
production  costs. 

Coal  is  transported  at  level  cT^  from  coal  region  i  to  demand 
region  j  .  Crude  oil  is  transported  from  oil  region  i  to  refinery  j 
in  a  quantity  oT^,  where  it  is  refined  into  a  fixed  proportion  of 


light  and  heavy  oil.  Refinery  1  converts  60%  of  its  crude  oil  into 
light  oil,  while  refinery  2  converts  50%  of  its  crude  oil  into  light 


oil.  Light  oil  L . .  and  heavy  oil  H.. 

LI  ij 

i  to  demand  region  j  . 


are  transported  from  refinery 


The  linear  program  representing  the  minimum  cost  of  production 
and  delivery  of  consumable  goods  to  meet  a  demand  q  e  is  of  the 
form 

min  <r,x>  subject  to 

Ax  =  q  demand  requirements 

Bx  =  0  transportation  flow  balance 

Ex  =  a  factor  of  production  upper  bounds 

Fx  =  b  steel  and  capital  resource  constraint 

x  =  0  non-negitive  factors  of  production 


The  components  of  x  represent  the  levels  of  production  of  coal  and 
crude  oil,  levels  of  coal  transported  to  demand  regions,  levels  of 
crude  oil  transported  to  refineries,  and  levels  of  refined  oil  trans¬ 
ported  to  the  demand  regions.  The  coefficient  r^  is  the  unit  cost 
of  the  activity  associated  with  x^.  The  vector  Ax  consists  of  levels 
of  consumable  goods  delivered  to  the  demand  regions.  Bx  represents 
the  conservation  of  materials  transported  through  a  region  or  re¬ 
finery.  Fx  is  the  vector  of  steel  and  capital  consumed  in  the 
processing  of  factors  x. 


The  demand  relationship  between  the  desired  quantity  of  a  con¬ 


sumable  good  in  a  particular  demand  region,  q^,  and  the  prices  con¬ 
sumers  are  willing  to  pay  for  such  goods,  { p ^  }  j=l,...,6,  is  given 


4 .  NEWTON'S  METHOD  AND  LEMKE'S  ALGORITHM 

An  equilibrium  problem  which  is  equivalent  to  that  of  Hogan's 

example  when  consumable  good  prices  and  quantities  are  positive  is 

the  following.  The  supplier  solves  the  linear  program 

A  >  *  >  > 

min  <c,x>  ,  subject  to  Ax=q  ,  Bx=b,  and  x=0  . 

Note  that  the  original  form  of  Hogan's  example  contains  transportation 
flow  equations  Bx=0.  However,  due  to  the  flow  conservation  nature 
of  the  constraints,  the  matrix  B  factors  into  the  form  B  =  (I,B). 

Hence,  replacing  0=Bx  with  -Bx2  >=  0,  where  x  =  (x^x^  and  x^=  -Bx2, 
and  similarly  decomposing  the  other  inequalities  and  replacing  x^ 
with  -Bx2>  one  obtains  a  linear  program  with  no  equality  constraints 
and  fewer  variables.  The  matrix  B  represents  all  non-demand  con¬ 
straints,  while  A  represents  the  demand  constraints.  The  subscript 
2  on  x  has  been  dropped.  The  supply  equals  demand  equilibrium  condi¬ 
tion  is  replaced  by  the  following  conditions: 

P-PD(q)  =  0  ,  q  =  P  ,  <q,p-pD(q)>  =  0  . 

These  conditions  are  that  no  supply  price  p^  is  ever  below  the  cor¬ 
responding  demand  price  pD(q)^,  and  if  a  positive  quantity  of  q^  is 
demanded,  then  p^  =  p^(q)^.  The  equilibrium  conditions  for  this  problem 
are,  dropping  hats  on  A,  B  and  c, 
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(1) 


T  T  >  >  T  T 

c-Ap-Bs=0,x  =  0,  <x,c-A  p  -  B  s>  =  0  ; 

(2)  Ax  -  q  =  0  ,  p  =  0  ,  <p,Ax  -  q>  =  0  ; 

(3)  Bx  -  b  =  0  ,  s  =  0  ,  <s,Bx-b>  =  0  ; 

(4)  P  -  PD(q)  =  0  ,  q  =  0  ,  <q,p-pD(q)>  =  0  . 


We  can  now  write  equilibrium  conditions  (1)  -  (4)  as  a  generalized 


F(z):  = 


e  Josephy 

[3]) 

.  Let 

K: 

II 

+  ^ 

where  d 

:  =  itd 

n 

n 

,  s  e  R 
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q  e  R 
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where  z:  =  (x,p,s,q).  Then  finding  an  economic  equilibrium  is  equiv¬ 
alent  to  finding  a  solution  to  the  generalized  equation 

(5)  0  e  F (z)  +  NR(z)  • 

The  linearization  of  (5)  at  a  point  z:  =  (x,p,s,q)  is  given  by 


(6) 

where 


LF-(z)  = 


0  e  LF-(z)  +  N  (z)  , 

2  K 


c 

r 

0 

-at 

T  1 

-B  0 

o 

A 

0 

0  -I 

-b 

♦ 

B 

0 

0  0 

-pD(q)D^(q)q 

0 

I 

o  -Pp(q) 

I 
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The  linear  complementarity  problem  corresponding  to  the  linearization 
(6)  about  z  is 

LF-(z)  =  0  ,  z  =  0  ,  <z,LF-(z)>  =  0  . 

z  z 

L eniK e  s  algorithm  (see  Cottle  and  Dantzig  [1])  can  be  aoplied  to  this 
linear  comnlementarity  problem  if  the  matrix 


— 

0 

-at 

-bt 

0 

A 

0 

0 

-I 

M:  = 

B 

0 

0 

0 

0 

I 

0 

-P^Cq)  _ 

is  copositive  plus,  where  copositive  plus  is  defined  as  follows: 

Definition  1.  A  matrix  M  is  copositive  plus  if  and  only  if 

(1)  <u,Mu>  =  0  for  all  u  =  0  ,  and 

(2)  if  u  =  0  and  <u,Mu>  =  0  ,  then  (M+M^)u  =  0 

See  Cottle  and  Dantzig  (  1  )  for  further  details.  Since  M  is  co¬ 

positive  plus  if  and  only  if  -p^Cq)  is  copositive  plus,  the  applic¬ 
ability  of  Lemke's  algorithm  rests  on  showing  that  -p^(q)  is  copositive 
plus.  Recall  the  demand  relations 

^n(qi/q°)  =  E  e..  In  (p./p?)  ,  i  =  l,...,6  . 

j  J  J  J 

Inverting  these  relations  yield 


where  (a)  =  (e) 


The  elasticity  matrix  c  is  of  the  form 


where  E:  =  (e^)  are  the  elasticities  of  quantity  versus  price  for 
coal,  light  oil  and  heavy  oil.  The  zero  elements  in  e  correspond 
to  the  independence  of  quantities  demanded  in  one  region  and  prices 

in  the  other  region.  If  we  define  F  =  (f„)  =  E  1  ,  then  it  is  easily 
seen  that 


The  matrix  -E  has  the  following  properties  for  Hogan's  example: 

(1)  positive  diagonals  , 

(2)  negative  off-diagonals  , 

(3)  strict  diagonal  dominance  , 


where  (Plemmons,  [5]  )  a  matrix  M  is  strictly  diagonally  dominant  if 


l 


and  only  if  there  exist  positive  constants  {d.}  such  that  for  each  i  , 


m..d.>  E  I  m .  .  ]  •  d  . 
11  1  jti  13  3 


For  the  Hogan  example,  the  d^^  can  be  chosen  to  be  1.  It  follows  from 


Plenunons  (  5  )  that  -F  is  strictly  positive.  The  matrix  -p^(q)  is 


(-PD^)ij  *  ^  “ij 

qj 


Thus,  for  positive  p  and  q,  -p^(q)  is  non-negative,  and  hence  Condition 


1  of  the  definition  of  a  copositive  plus  matrix  is  satisfied.  To 
show  that  Condition  2  also  holds,  it  will  be  convenient  to  define 


-pi 


a  matrix  G  =  -p^(q),  where  G. .  =  ~  a.  .  .  Thus,  G  has  the  form 


qj 
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v/herc  all  g  >  0  ,  i=l,...,6  ,  j=l,...,3 


To  show  the  validity  of  Condition  2  of  Definition  1,  suppose 
u  =  0  and  <u,Gu>  =  0  .  Let  w:  =  (u  ,u_,u  )  and  z:  =  (u  ,u.,u  ). 


-i0- 


.  I 
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-r6— - 
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Let  gi  denote  the  row  fgil  '2i2'Si3) ,  and  g^-w:  = 

RilKl+8i2W2+gi3U3  ’  for  i=1»---»6-  Thcn  Gu  =  (gj . 'v>&2' ‘Z’g3. 'W' 
g4*'2’  g5*  W,g6*  and  <u>Gu>  ~  U1*(S].*W)  *  u2*(g2#*z)  + 

U3*  (g3*  +U4  *  (g4  •  ’  Z^+U5*  *w)+u£)*  ("(. .  * z)  •  Suppose  some  u^O, 

say  Uj.  Then  >  0  and  thus  g^u^  >  0  ;  hence  g  *w  >  0  and  thus 
ul*(8i.*'0  >  0,  contradicting  0  =  <u,Gu>.  A  similar  argument  works 
for  any  choice  of  >  0.  Hence,  if  0  -  <u,Gu>  ,  then  u  =  0  .  Thus, 
Condition  2  holds,  and  -p^Cq)  is  copositive  plus.  Indeed,  the  above 
analysis  actually  shows  that  -Pp(q)  is  strictly  copositive.  (Sec 
Cottle  and  Dantzig  (  1  J  for  details  on  strict  copositivity). 

An  interesting  property  of  -p^fq)  relates  strict  diagonal 
dominance  of  -p^(q)  to  the  same  property  for  -a  .  Note  that 


■pD  (q)  ij 


"£i  0 

qj  ij 


For  all  q^  >  0,  take  d_.  =  q^  in  the  definition  given  for  strict  diagonal 
dominance.  Then  -Pp(q)  is  strictly  diagonally  dominant  if  and  only 
if,  for  each  i, 


For  p^  >  0,  this  reduces  to 


«  !i 

q j  ii  i 

.i^i  qj 

d .  with  q.  , 

we  have 

l  *i 

-piaii  ”  pi 

•  E  |a 

reduces  to 

ij 


-  a. .  >  E  a . . 
11  j*i  13 
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Thus,  -p^(q)  will  be  strictly  diagonally  dominant  when  -a  is  strictly 
diagonally  dominant.  The  implications  of  strict  diagonal  dominance 
are  discussed  in  great  detail  in  Plemmons  (  5  ). 

A  final  remark  is  in  order  concerning  the  assumption  of  posi¬ 
tive  prices  and  quantities  of  consumable  goods.  The  form  of  the 
function  p^  assumed  by  Hogan  is  undefined  when  any  quantity  q^  is 
zero,  since  the  diagonal  elements  of  a  are  negative.  Thus,  an  equi¬ 
librium  solution  which  satisfies  the  condition  that  the  supplier's 
price  must  not  be  less  than  the  consumer's  price  must  have  q  >  0 
and  (since  PD(q)  >  0)  p  =  PD(q)  >  0.  Thus,  it  follows  from  Josephy 
[3,  Theorem  1]  that  in  a  sufficiently  small  neighborhood  of  the  equi¬ 
librium  nrice  and  quantity,  the  linearized  problem  will  yield  positive 
p  and  q  . 
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Newton's  method  for  generalized  equations  (Josephy  [3])  has  been  applied  to 
the  economic  equilibrium  problem  of  the  Proiect  Independence  Evaluation  System 
(PIES)  Energy  Model  (Josephy  [4]).  The  resulting  algorithm  involves  solving  a 
sequence  of  linear  complementarity  problems.  Lemke's  complementary  pivot  algo¬ 
rithm  is  used  for  this  purpose.  In  this  paper,  we  show  that  the  linear  comple¬ 
mentarity  problems  will  be  copositive  plus  when  the  negative  of  the  elasticity 
matrix,  -e,  of  the  consumer's  quantity  vs.  price  relation  has  the  following  prop 
erties:  (i)  positive  diagonals,  (ii)  negative  off-diagonals,  and  (iii)  strict 
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Abstract  (continued) 


diagonal  dominance.  These  conditions  are  satisfied  for  Hogan's  example 
(ilogan  [ 2  ] ) .  Thus,  Lemke's  algorithm  will  either  converge  to  a  solution 
or  show  that  no  solution  exists.  Under  the  conditions  of  Josephy  [3,  Theorem  1], 
a  solution  to  the  linear  complementarity  problems  will  always  exist.  Hence, 
lemke's  algorithm  can  be  used  when  the  conditions  of  the  Theorem  1  of  Josephy 
I  3 1  are  satisfied. 


